Bulk dynamics for interfacial growth models 
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We study the influence of the bulk dynamics of a growing cluster of particles on the properties 
of its interface. First, we define a general bulk growth model by means of a continuum Master 
equation for the evolution of the bulk density field. This general model just considers arbitrary 
addition of particles (though it can be easily generalized to consider substraction) with no other 
physical restriction. The corresponding Langevin equation for this bulk density field is derived where 
the influence of the bulk dynamics is explicitly shown. Finally, when it is assumed a well-defined 
interface for the growing cluster, the Langevin equation for the height field of this interface for some 
particular bulk dynamics is written. In particular, we obtain the celebrated Kardar-Parisi-Zhang 
(KPZ) equation. A Monte Carlo simulation illustrates the theoretical results. 
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I. INTRODUCTION 

In the last fifteen years there has been a great inter- 
est in the study of the growth of surfaces by dynamic 
processes based in addition and substraction of particles 
(see, for example, Jl|-Ql). For instance, the understanding 
of the conditions under which a growing surface shows a 
rough structure is nowadays of the greatest importance 
in the production of thin films and/or pure cristals. Sur- 
face growth is usually studied by using lattice models in 
which simple stochastic rules intend to mimic the rele- 
vant phenomena. Their extensive computer simulation 
have been of a major importance in characterizing and 
understanding the different morphologies that occur in 
real experiments. However, due to the intrinsic limi- 
tations of computers capabilities, some very interesting 
aspects are usually subject to inconclusive analysis and 
data interpretation. In particular, let us remark the in- 
herent difficulty in the study of the surface long time be- 
havior and its scaling properties. Nevertheless, for this 
particular aspect, analytical models seem to give us the 
answers to the questions that the computers fail to clar- 
ify. These are mainly based in postulating a Langevin 
equation for the height of the interface measured from a 
reference substrate. Such Langevin equations intend to 
mimic the system microscopic dynamics and its collective 
effect at large scales in space and time. A general choice 
has the structure 



+ ... + %(*) , 



(1.1) 



where /it(x) is the height of the interface at time t at 
the substrate position x G R d and 7?t(x) is a white noise 
term. It is generally assumed a on e-to-one correspon- 
dence between various terms in (1.1) and different phys- 
ical processes (for example see the discusion in || con- 
cerning a model for epitaxial growth, or a general method 



in H to propose an equation such as (1.1) by using the 
reparametrization invariance symmetry). The details of 
the microscopic processes that are assumed to be irrele- 
vant at this scale of observation are taken into account 
through the values of the coefficients v\ , 1/2, ■ ■ ■ and 
the prope rties of the noise term r\ t ■ Once the Langevin 
equation (1.1) is defined, one may apply renormaliza- 
tion group procedures to obtain different universal prop- 
erties and scaling behaviors. The success of this scheme 
is clearly represented by the definition and analysis of 
the Kardar-Parisi-Zhang equation (KPZ) which has 
been a clear breakthrough in the study of the space-time 
asymptotic behavior of growth models. 

Quite often surface growth is a consequence of bulk dy- 
namic processes. A good example of this is provided by 
the growth of bacterial colonies where bacteria multiply 
in a nutrient environment, the shape of the colony being 
the moving interface || . In general, the dynamics of the 
particles before and after its aggregation to a substrate 
may influence the system interfacial behavior. For in- 
stance, it may appear shadowing effects as it happens in 
DLA processes and in thin film growth j|, or they may 
induce different scaling regimes depending on the time 
interval studied as it is shown in some molecular beam 
epitaxy models in which a system bulk dynamics is de- 
fined [po|JTT] . However, interface models are usually ex- 
pressed in terms of a height field, h t (x) > 0. In doing so, 
bulklike contributions are neglected since only interfacial 
degrees of freedom are being considered. Unfortunately, 
the mathematical hurdles to deriving the phcnomcnolog- 
ical dynamics of interfaces from stochastic bulk micro- 
scopic models is formidable, and a comprehensive theo- 
retical picture is still lacking although a significant body 
of rather rigorous work has been devoted to the subject 
& 

It is well known that the macroscopic behavior of sys- 
tems at nonequilibrium states exhibits a strong depen- 
dence on the functional structure of its microscopic dy- 
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namics (for instance, in the so-called two-temperature 
Ising model |l3|] one finds that the phase diagram changes 
radically depending on the analytic form for the proba- 
bility of a spin- flip). Nevertheless, one may expect that 
this strong relation between microscopic dynamics and 
macroscopic behavior should dissappear near a renormal- 
ization group (RG) fixed point or in the scaling regime 
where universality seems to guarantee that the micro- 
scopic details are irrelevant (at least, one knows that this 
is true when studying dynamic properties of equilibrium 
systems near a (RG) fixed point However, some 

recent results on the critical behavior of a nonequilib- 
rium driven diffusive system show that the microscopic 
dynamics may play a relevant role in the determination 
of the system universality class . The influence of the 
microsocopic dynamical details into the critical and non- 
critical properties of a nonequilibrium model implies, in 
our opinion, that any a priori construction of a Langevin 
equation as (l.f ) may occasionally disregard important 
features. 

In this paper we introduce a quite general class of 
nonequilibrium bulk growth models for which the afore- 
mentioned problems can be addressed (there are in the 
literature some efforts in this direction |l6|-|l8|] ) . We shall 
define a stochastic bulk local dynamics expressed by an 
appropiated continuum Master equation in which, for 
simplicity, only addition of particles is considered. From 
the Master equation and using a truncated Kramers- 
Moyal expansion, we shall derive a Langevin equation for 
the bulk degrees of freedom in which there is an explicit 
dependence on the analytic form of the rates. In order to 
study surface properties, in a subsequent section we shall 
derive, from this bulk equation, an expression for the in- 
terfacial height field dynamics. The illustrative example 
we take is that of the KPZ equation and, for consistency, 
in that particular case we check our results by means of a 
Monte Carlo simulation. Some other examples are then 
briefly commented and our conclusions are given in the 
final section. 



II. GROWTH MODELS: A GENERAL 
DESCRIPTION 

Let us consider a particle density field, p(x,r), x 6 
R d+1 , and assume that the probability distribution P T {p) 
associated with each field configuration obeys the follow- 
ing continuous Master equation 



d T Pr{p) = / dv 
l R d+l 



c( Pr ->■ p)P T {Pv) 

-c(p - p T )Pr(p) 



(2.1) 



Here, c(p — > p') is the probability per unit time (or tran- 
sition rate) from one configuration p to another p' , and 
p r = p(x) - fT^x - r), p r = p(x) + Q-^ix - r). Note 
that the density field, p, can only grow in steps of size 



Q , This is consistent with a picture in which p(r) is 
a particle density that results after coarse graining over 
blocks of size SI centered around r in a lattice. Therefore, 
the Master equation (2.1), so defined, could be thought 
of as if it only described processes that add one particle 
per block of the lattice per elementary time step. This is 
schematically represented in figure [j]. 

Next, we choose the transition rates such that 



c(p -> p ) = w(p;r), 



(2.2) 



namely, they are a function that depends only on the 
initial configuration p and the specific point where mass 
is added. Now, let us assume that is large enough 
(tipically, it should be much bigger than any microcopic 
length scale present in the original physical problem but 
much smaller than the correlation lcn ght of the system) 
and expand the Master equation (2.1) on invers powers 
of f2. Then, using the expansions 



C (P -> P r ) = w(p;r), 

c(pr —* p) = w[p — fl~ S(x — r); r) 
(-l) m S m 



E 



m=0 PV ; 



(-1) 



m\tt m Sp{r) 



w(p; r), 



-PAP), 



(2.3) 



(2.4) 



(2.5) 



we get the so-ca lled Kramers-Moyal expansion of the 
Master equation (|2~T|) [111, 



drPAp) 



-i) 



l si 



^ HQ 1 Sp(r) 1 



w(p;r)P T {p) 



(2.6) 



The ne xt s tep we take is to keep only the first two 
terms in (2.6). Then, we can write down the following 



Fokker-Planck equation 



d T P T {p) = f dr 

JB d+1 



1 



SI Sp(r) 

1 s 2 



2tl 2 5p(r) 2 



(w(p;r)P t (p)). (2.7) 



To control the goodness of such approximation we appeal 
to Kurtz theorem pc| , by virtue of which when fl — > oo, 
and for a given time T < oo then 



_ >r logf2 

SUp \p T - p T \ < Csv 



V. 



(2.8) 



where p T and p T are typical time trajectories on phase 
spac e which are solutions of the exa ct Master equation 
( |2.lD and the Fokker-Planck one (2/7) respectively. is 
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a random variable whose distribution does not depend on 
f2 and satisfying (exp(A^Q )) < oo for any constant A > 0. 
That is, for a given fixed time T one can always find a 
large enough fl such that the difference between solv- 
ing exactly the master equation or solving the truncated 
version of it, is of order logf2/f2 and therefore negligible. 
Moreover, this bound is the best one and no new terms of 
the Kramers-Moyal expansion give better results. How- 
ever, when T — > oo one cannot control, in general, the 
accumulated influence of the neglected terms in the ex- 
pansion. But, since the study of growth models is manily 
focused on the undestanding of their evolution properties, 
due to Kurtz's theorem the Fokker-Planck equation (2.7) 
is a valid theoretical starting point. 

Lastly, the Fokker-Planck equation ( |2.7| ) is equivalent 
in the Stratonovich sense to the Langevin equation p9| 



p(r,t) = Q(h t (y 



(3.1) 



d t p{r,t) = w(p;r) - 
1 



4ft 6p(r) 
w(p;r)V 2 v t (r) 



w{p;r) 



(2.9) 



where t = fl r, Vt(r) is a white noise with zero mean 
and {v t {v)v t ,{r')) = S(r - r')8(t - t'). We recall that 
w(p; r) is the probability per unit time of adding a par- 
ticle (of mass f2 _1 ) at point r. If we considered the posi- 
bility of particle substractions the n we should do the fol- 
lowing substitutions in equation (2.9): w — > w+ — tw- 
in the first term and and w — > w + + u>_ in the last two 
terms. w+(-) = w(p', r )+(-) i s the probability per unit 
time of adding (substracting) a particle. 

Before we pro ceed to the next section, let us remark 
that a) eq. (2J)) describes the evolution of the bulk den- 
sity field of a system that growths by addition of particles 
with, in principle, no other physical assumptions, and 
b) the Langevin equation (2.9) depends directly on the 
functional form of the bulk rates. The election of these 
bulk rates then provides the physical restrictions for the 
particular growth model that is going to be specifically 
modelled. Also, it is remarkable that the influence of the 
bulk dynamics affects the noise term through a nontrivial 
factor. 



III. HEIGHT DYNAMICS: KPZ EQUATION AND 
MONTE CARLO SIMULATION. 

We proceed to single out the interfacial degrees of free- 
dom of equation (2.E). In order to achieve th is, w e place 
two conditions on the solutions of equation (2J5). First, 
let us impose that our bulk dynamics produces a surface 
perpendicular to the z-axis without overhangs and bub- 
bles. This condition is necessary to ensure t hat we have a 
well defined interface (note that equation (|2.9|) contains 
overhang/vacancy and shadowing effects). Secondly, we 
shall neglect any interfacial profile. Then, we may as- 
sume that the solutions of the Langevin equation ( |2.9| ) 
have the form 



where r = (x, z), x 6 R d is a point in the substrate 
and /i*(x) is the height of the growing surface at time 
t. 6(A) = 1, 1/2, if A > 0, = 0, < respectively. That 
is, for a given point in the bulk r, if its z coordinate is 
larger, equal or smaller than the actual position of the 
surface, h t (x.), then the density field is p(r,t) = 0, 1/2 
and 1 respectively. Note that since the function is not 
continuous, when differentiating we should use a regu- 
larized version of it, e.g. Q(x) = + tanh(aa;)] with 
a — » oo. 

We are interested in constructing a dynamical equa- 
tion for the h t (x) field s. T herefore, let us make a time 
derivative in equation (|3.l|) 



d t p{r,t) = d t ht(x)5(ht(x) - z). 



(3.2) 



Integrating in z both sides of equation (3.2) we find 



d t h t 



dz d t p(r,t). 



(3.3) 



Equating this expression for d t h t (x), together with (pTs|) , 
will lead us to the desired Langevin equation for the 
heights. 

To make this a bit more concrete, we now introduce a 
particular set of rates. For instance, we choose the prob- 
ability of adding mass to the point r to be proportional 
to the square of the gradient of the density field in that 
point: w(p;r) oc |Vp| 2 . With this election the unwanted 
effect of nucleation of bubbles is avoided. After a bit of 
algebra, we get 



d t ht 



a(l + (V/i f ) 2 ) + 



D 
2ft 



Ah 



ft 1 / 2 



l + (Vh t Y) ut, 



1/2 



(3.4) 



which is the celebrated KPZ equation with a different 
noise term (a naive po wer counting argument reduces 
the relevant part of (3.4) to the KPZ equation). The co- 
efficient D has the proper dimensions and a is positive 
and depends on how the O function is regularized. This 
comes from our particular ansatz, but, we would like to 
stress here that, as far as universal properties are con- 
cerned, the precise value of the coefficients is immaterial. 
In fact, it is easy to show (with naive power counting) 
that for any bulk dynamics given by w(p;r) oc |V j o|' ? 
with rj > 0, gives rise to a height equation falling in the 
KPZ universality class. 

Next, we proc eed to check numerically the connection 
betw een eq.( |2.l|) with w(p; r) oc | Vp| 2 and the KPZ equa- 
tion fljTg). 

Numerical results. The simple bulk rate |Vp| 2 can be 
easily implemented in a Monte Carlo experiment. On a 
two-dimensional square lattice periodic boundary condi- 
tions are considered in one of the principal axes. Each 
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lattice site is labeled by an occupation variable p r rang- 
ing from 0, l/f2, ... to 1. A site is empty if p r = and 
full if p r = 1. Initially the lattice is empty except for a 
full horizontal bottom line. The growth starts when an 
empty site r is chosen at random from the lattice. Then, 
p r is increased in fl^ 1 with a probability that depends 
on each nearest neighbour of r through |Vp r | 2 . Since 
symmetrized discrete forms of the gradient operator may 
give rise to bulk vacancies incompatible with ( |3.l[ ), we 
use the following finite-difference formula for the density 
derivatives (r = (x, z)) 



-c(p -> ff)P T {p) 



(^p{x ± 1, z) — p(x, z), p(x, z) — p(x, z — 1) j . 



(3.5) 



That is, left and right derivatives are used alternately 
to avoid asymmetric effects and, for convenience, a unit 
lattice spacing is assumed. The height is defined as the 
distance to the highest occupied lattice site directly above 
the substrate coordinate x. 

Figure |2| shows the scaling plot for the surface width 
W(L, t) 2 = L- l JJ]H{x)-h{t)] 2 for different system sub- 
strate sizes L. h(t) is the mean height of the interface at 
time t. The numerical data were averaged over 2000 in- 
dependent runs for L — 100, 200, 500 and over 1000 inde- 
pendent runs for L — 800. For reasons of computational 
efficiency the results shown correnspond to CI = 1 . Other 
values of fl yield similar results but the simulations are 
much more time-consuming. Good data colapse is ob- 
tained for a roughness exponent a = 1/2 and a dynamic 
exponent z — 3/2, in agreement with the KPZ prediction 

0- 

The example we have just provided is by no means 
unique. Our formalism emcompasses many others well 
known growth equations. Let us just mention that with 
the simple dynamics given by w(p,r) = |Vp| the equa- 
tion of Golubovic and Wang, related to the anharmonic 
equilibrium thermal fluctuations of smectics A [pl| , is 
obtained. This is given by 



d t ht 



1 



+ 



i + (y ht y)(\ + -H) 

1/4 

l + (X7h t ) 2 ) v t , 



(3.6) 
(3.7) 



where A is a coefficient and H is the curvature (see [^l| ) . 
As we have mentioned before, this kind of dynamics (pro- 
portional to \V pY 1 with r] > 0), falls in the KPZ univer- 
sality class. Also, the Edwards- Wilkinson equation [^2| . 
which favours growth at local minima, can also be recov- 
ered by considering substraction of particles and a rate 
of the form w — |V 2 p|. In this case, the formalism has 
to be slightly modified by linking the election of additon 
or substraction of particles to the density field configura- 
tion. More explicitely, now the Master equation defining 
the process reads 



drPr(p) 



dr 



(3.8) 



with p' = p(x) + aJ7 _1 <5(x — r) and a — —1,0,1 for 
V 2 p less, equal and greater than 0, respectively. That is, 
material is added to those areas where the laplacian of 
the density field is negative and taken from those where 
it is positive. In this manner, a balanced distribution of 
mass is achieved that results in the equilibrium Edwards- 
Wilkinson universality class. 

Many other different rates lead to their corresponding 
growth equations, sometimes to the same one, showing 
that growth models with similar surface behavior may 
not have the same bulk properties. 



IV. CONCLUSIONS 

In this paper we have introduced a class of nonequi- 
librium models in which a stochastic bulk dynamics is 
defined. The bulk evolves by an adsorption process rep- 
resented by a continuum Master equation. From it, we 
have derived a Langevin equation for the bulk density 
field whose structure depends on the details of the un- 
derlying bulk dynamics. This dependence was then ex- 
tended to an equation of motion for the interfacial degrees 
of freedom. In particular, we have examplified the pro- 
cedure by deriving the KPZ equation from a very simple 
bulk rate. A Monte Carlo simulation of the bulk process 
confirms the predicted scaling behavior for the interface. 
Finally, a number of examples were briefly mentioned. In 
all cases the bulk dynamics determines the mesoscopic 
height equation, showing that both scales are related in 
a non trivial form and that their mutual influence could 
be far from intuitive. 

The strategy developed in this paper is quite gen- 
eral. It includes both local and nonlocal, and equilibrium 
and nonequilibrium growth processes. Therefore, a great 
number of growth physical phenomena can be studied, 
in principle, with our approach. For instance, molecular 
beam epitaxy models with adatom mobility, driven lat- 
tice gases or wetting phenomena by means of lattice gas 
theories of multilayer adsorption, to name just a few. 
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FIG. 1. Schematic rep rese ntation of the proccess described 
by the Master equation (2.1). 
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FIG. 2. Scaling plot for the model defined by the Master 
equation (2.1) with the rate |Vp| 2 . The data ara for L = 100, 
200, 500 and 800, where L is the width of ths substrate. 
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